% 创建时间：2025-06-20 22:55:50
% 作者：兪馳
clc;
clear;
A=[1,3,5,7;
   2,4,6,8];
[L,A1,M]=svd(A),
A2=L*A1*M',
norm(A-A2)